Detecting phase synchronization by locahzed maps: 
Apphcation to neural networks 
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We present an approach which enables to state about the existence of phase synchronization in 
coupled chaotic oscillators without having to measure the phase. This is done by observing the 
oscillators at special times, and analyzing whether this set of points is localized. In particular, we 
show that this approach is fruitful to analyze the onset of phase synchronization in chaotic attractors 
whose phases are not well defined, as well as, in networks of non-identical spiking/bursting neurons 
connected by chemical synapses. 



Many neural networks rely on a synchronous behavior 
for a proper functioning, e.g. information transmission 
P, Q , pattern recognition 3] , and learning 4| . Never- 
theless, the onset of synchronization in a network may 
also lead to some diseases as Parkinson disease Q and 
epilepsy Q . Studies on this topic have been concentrated 
on synchronization of networks of identical chaotic oscil- 
lators, in which the onset of complete synchronization 
takes place However, they do not treat the onset 
of synchronization in networks of non-identical chaotic 
oscillators, case typically found in nature 8], where com- 
plete synchronization is much harder to be achieved. In- 
deed, in such networks a weaker kind of synchronization 
may take place, the phase synchronization (PS). The con- 
dition for PS between two subsystems k and j can be 
written as: 



\(l)k{t) - r(j)j{t)\ < c, 



(1) 



where (j)k,j is the phase of the subsystem Sfcj-, c S M is 
a constant, and r is a rational number. PS is a com- 
mon phenomenon in interacting chaotic oscillators 
and plays a major role in physical processes linked to 



communication 
human brain 

In order to state about the existence of PS, one has 
to introduce a phase 4>(t) for the chaotic oscillator, what 
is not straightforward. Indeed, there is no general defi- 
nition of phase to chaotic attractors, and depending on 
the situation one has to decide which phase is more suit- 
able. In fact, in some oscillators it is rather unclear which 
phase one should use, especially in non-coherent oscilla- 
tors with more than one time scale, typically found in 
neuronal dynamics with bursting/spiking behavior. 

In this letter, we present a general and easy way to 
identify PS without having to access explicitly the phase. 
The approach consists in defining maps, which are a 
natural extension of the stroboscopic map, to coupled 
chaotic oscillators, in which the oscillators are observed 
at especial times. PS implies the existence of maps of the 
attractor that appear as localized structures in the acces- 
sible phase-space. The fact that PS produces subsets of 
the attractor that are localized structures, by particular 



observations was previously used as a way to detect PS 
in chaotic oscillators [1, [111, [l2|- Here, we extent these 
results by demonstrating that localized sets can be con- 
structed while in PS by means of any typical physical 
observation, which has a strong impact in the field of ex- 
perimental physics, since in the laboratory measurements 
are restricted to the limitations of the experiment. Note 
that since this approach does not requite any further cal- 
culation, but just the analysis of whether the sets are 
localized, it can be used in real time experiments for PS 
detection. We illustrate the power of this approach by 
analyzing PS in a network of non-identical Hindmarsh- 
Rose (HR) spiking/bursting neurons connected via chem- 
ical synapses. 

The classic stroboscopic map is defined in periodically 
driven chaotic oscillators. It consists in sampling the 
chaotic trajectory at times hTq, where n is an integer 
and To is the period of the driver. The stroboscopic map 
was used to detect PS The basic idea is that if 

the stroboscopic map is localized in the attractor, PS is 
present. To generalize the stroboscopic map to coupled 
chaotic oscillators, we do the following: Given two sub- 
systems Efe and Sj, we observe Sfc at times when some 



^and communication processes in the event in the oscillator happens. As a consequence 

of these observations, we get a discrete set Vk- Then, 
we demonstrate that if there is PS then the set Vk is 
localized. 

In order to introduce our ideas, we analyze PS in two 
coupled chaotic oscillators, namely the Lorenz oscillator 
driven by the Rossler. Further, we extent this result to 
general compact oscillators. The subsystem corre- 
sponds to the Rossler oscillator and to the Lorenz. 
They are coupled unidirectionally in the driver response 
scheme. An event in the is considered to happen when 
its trajectory crosses a Poincare section yk = 0. As a re- 
sult, we get the series of times {rDi^fi, where t^. is the 
time at which the zth crossing of the trajectory of E^ 
occurs in a Poincare plane. The two coupled oscillators 
are given by: 



ifc = -a(j/fc + Zk) 
yk = a{xk + 0-2yk] 
Zk = a[0.2 + Zk{xk 

with 



5.7)] 

= 13,cr 16, 



^3 ~ ^iVi 



Vj ■''J ^3 

^I3z, 



= 45.92, and /3 = 4. 



Since the trajectory of the Rossler oscillator rotates 
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around a fixed point [Fig. [TJa)], we can define a phase 
Ok = tan^^{yklxk) which provides: 6k{t) = J^iykXk - 
Xkyk)/{xf + vDdt. The trajectory of the Lorenz does 
not have an unique center of rotation, see Fig. [Ijb). 
However, if we consider the projection (u,Zj) with u = 

■\Jx'j + yj , the trajectory projected into this subspace 
presents an unique center of rotation. Tlius wc also 
define a phase 6j — tan~^[{zj — ZjQ)/{u — uq)], where 
(Mo,2jo) = (19,45) is the center of rotation in the sub- 
space {u,Zj), which provides 9j{t) — JqIzju — u{zk — 
■2jo)]/[(w - uof + {zj - Zjof]dt 

For e = 0.0, we construct the set Vj by sampling the 
trajectory of Ej at times r^. The set Vj spreads over 
the trajectory of E^; there is no PS, the phase difference 
= 9k — 9j diverges [fig. [U (c,e)]. Indeed, a calcu- 
lation of the frequencies shows that {9k) ~ 13.94 and 
{9j) « 13.75. As we increase the coupling, PS appears. 
In particular, for e = 13.0 the set Vj is localized, and 
the phase difference is bounded [fig. [Hd,f)]. The average 
frequencies are {9k) = {9j) ~ 13.95. 





FIG. 1: PS onset between a Lorenz oscillator driven by a 
Rossler oscillator. The attractor of the oscillators are depicted 
in gray (a-d). In (a) we show the projection of the Rossler 
attractor into the subspace {xk, Vk)- The black line shows the 
Poincare section ai yk = with the constrain yk < 0. The 
projection of the Lorenz attractor in the subspace {xj,Zj) is 
depicted in (b), and in (c,d) into the subspace {u,Zj). For 
e — 0.0, the set Vj (in black) spreads over the attractor of 
the oscillator Ej (c), and there is no PS; the phase difference 
diverges (e). For e = 13.0 the oscillators present PS. The set 
Dj is localized, it does not fulfill the attractor of Ej (d); the 
phase difference is bounded (f). 

Next, we demonstrate that the sets Vj of the attractor 
that appear as localized structures imply PS, and vice- 
versa. We first show for a Poincare section for a better 



understanding of the ideas, and then we generalize these 
results to any possible event. S-,- is given by the dy- 
namical system Xj ~ Gj(xj), let Fj be the fiow, Tj the 
Poincare section, and IIj the Poincare map associated 
to the section Tj, such that given a point x^ S Tj, so 

x}+i = nj(x}) = F^"'''^\x)), where At^ ^ - t^~\ 
From now on, wc use a rescaled time t' — t/{Tj). For a 
slight abuse of notation we omit the /. The average return 
time is given by {Tk) = E»=o ^^4/^ = 1^^ and the 
time is rescaled, such that (Tfe) = 1. From our hypothesis 
that both oscillators present phase-coherent oscillations, 
there is a number Kk such that [r^ — i{Tk)\ < k^, where 
Kk ^ 1. If both oscillators are in PS, then {Tk) ~ {Tj), 
and so: 

H-r}\<ii, (2) 

with K < Kk + Kj 1. Now, we analyze one typical os- 
cillation, using the basic concept of recurrence. Given 
the following starting points x^ € Tk and x° € Tj, 
we evolve both until x^ returns to Tj. Let us intro- 
duce At' = AtI - Arj. So, F'^'h^") = n,(xO) = 
x] e P,. Analogously, F^^^^l) = pf "^'+^"' (x^) = 

F^^Ff Mx2J = F^\nfe(xO)) = F^^(xi). Now, by 
using the fact that |At*| < k, we can write: F^'^ (xj.) « 
x^ + G(x^)k + 0{k?). So, given a point x^ G P^ eval- 
uated by the time when the trajectory of Ej returns to 
the section Tj, the point x^ returns near the section Ffc, 
and vice-versa. For a general case, we have to show that 
a point, in the section F^, evolved by the flow for an 
arbitrary number of events in the oscillator Sj, still re- 
mains close to Tk- But, this is straightforward, since 
Yl!i=Q ^■''N — I'T'k ~ < ^- demonstrated that 

the PS regime implies the localization of the set Vk- Now, 
we show that the localization of the set T>k implies PS. 
Supposing that we have a localized set T>k, so, eq. ^ is 
valid, by the above arguments. Therefore, we just have 
to show that eq. ([2]) implies PS. To do so, we note that at 
every crossing of the trajectory with the Poincare section 
the phase increases 27r, as a consequence 0a;(t^) = 2i x tt. 
Then, \Mrl) - <I>,{^1)\ = \Mrl) - ^r] -f C)|, where 
C, = T^ — Tj. Now, expanding the phase (jjj in Taylor series 

around r], we have (pjirj+Q « 2i7r-|-(/.j(r]) xC + 0{C^), 
as a result, the phase difference can be written as : 

l'/'feK)-0,W)| <Ax |T^-r;| <Ax^, (3) 

where, A — maxtj{(t>j{t)}. Therefore, we showed that 
boundness in eq. Q implies a bound in the phase 
difference at times r^. However, since the phase de- 
pends smoothly on time, and the Poincare section can be 
smoothly changed, the boundness in eq. ([3|) also holds 
at the continuous time. Thus, we conclude our result. 

An important point to stress is that it is not always 
possible to define a Poincare Section on the attractor in 
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such a way that a phase increases 27r every crossing. As 
an example we quote the non-coherent attractors with 
no proper rotations, where the definition of such section 
is not possible. Moreover, even if the oscillators are co- 
herent, it might happen that the accessible data is not 
suitable to define a section, but rather to measure the en- 
trance of the trajectory in some small region of the phase 
space. That does not constitute a problem, because PS 
implies localization of the set "Dfc, independently on the 
event definition. 

Let us first discuss the idea of localization. If the set T) 
is a subset of we say that V is localized (with respect 
to $) if there is a cross section ^ and a neighborhood A 
of ^, such that I? U A = 0. In particular, for practical 
detections, one may check whether V is localized, by the 
following technique. If there is PS, for y G "D it exists 
infinitely many x S A such that y fl -B^(x) = 0, where 
i3^(x) is an open ball of radius £ centered at the point x, 
and i is small. Then, we may vary y, x (one may take 
X to be an arbitrary point of the attractor) and £, to 
determine whether T> is localized. 

The event definition that generates the time series 
{TfcjiGN can be arbitrary. Therefore, the event could be 
a local maximum/minimum, the crossing of a dynamical 
variable with a threshold, the entrance in an e-ball, and 
so on. The only constraint is the event must be typi- 
cal. We also suppose that there is a function phase 0^, 
in such a way that (j)k = where Q,k is continuous 
and fifc < T. Under such hypotheses, we can state that: 
Given any typical event, with positive measure, in the os- 
cillator T,kj, generating the times (t^^ ) igN, if there is PS 
the observation of Yi^ at (tp^gN generates a localized set 

Next, we demonstrate this result, doing so we extent 
these ideas to non-coherent oscillators. The strategy to 
demonstrate the previous results to an arbitrary event is 
the following: (i) Note that the phase 4>j{t) naturally de- 
fines a section, namely Fj, in the attractor such that at 
the A^th crossing of the trajectory of Ej with Pj the phase 
is equal to A^ x 2ti. Obviously, this section depends on 
the initial conditions, (m) Suppose that we construct the 
set Vk by observing the trajectory of Efc at every cross- 
ing of the trajectory of Sj with Pj. Then, following the 
previous results, PS implies the localization of 2?^, and 
vice- versa. {Hi) Suppose that we have a small piece . 

of the section P^, such that the crossings of the trajectory 
of Yij with Pp produces a subsequence (T"')„;gN of the 
sequence (Tj)igN- Thus, we just note that if the observa- 
tion of the trajectory of at times (Tj)igN gives place 
to a localized set Vk, the observation at times (T"')„.gN 
also gives place to a localized set Vk which is a subset 
of Vk- Therefore, we showed that the observation of the 
trajectory of E^, when the trajectory of Ej returns to 
Pj=.,, also leads to a localized set in E^. {iv) Next, we 
show that an event does not have to be a piece of the 



section f ^ in order to obtain a localized set in E^. In- 
deed, given an e-ball event that produces the time series 
fj, in Ej, there is, at least, one intersection of this ball 

with the section Pj . Since Pj depends on the initial con- 
ditions, we can choose an initial condition right at the 
£-ball event. Next, we choose Pf such that it is com- 
pletely covered by the e-ball. Since the measure of the 
e-ball is small, e <C 1, the time difference between cross- 
ings of the trajectory with Pf . and the e-ball, namely 

fj — r"' , is also small. Therefore, if we observe the tra- 
jectory of Efc at times (fj)igN, we get a localized set in 
Efe close to the set T>k- Thus, we conclude our result. 

In order to illustrate these ideas, we consider again the 
Lorenz oscillator driven by the Rossler oscillator. As we 
showed before, for e — 13.0 there is PS. Thus the sets V 
must be localized independently on the event chosen. We 
define the event in the oscillator E^ to be the crossing of 
the trajectory 5fe = {xk,yk,Zk G R\xj < -8,yk = 0,yk < 
0}. These crossings generate the times (t^)^^^. Sk is 
depicted in black bold line in Fig. ^a.) together with the 
attractor of the Rossler oscillator depicted in gray. The 
observation of the trajectory of Ej at the times (tfe)ieN 
generates a localized set Pj[Fig. (Hb)]. 
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FIG. 2: PS implies the existence of localized sets, which is 
generated by the observation of an arbitrary typical event. 
The attractors of the oscillators are depicted in gray. For e — 
13.0, we show the projections of the Rossler Lorenz attractors 
(a,b). The black bold line shows the segment Sk (a), ©j 
is constructed by observing the oscillator Ej whenever the 
trajectory of Efe crosses the segment Sk- In (b) we show the 
set Dj in black. Since there is PS, fj is localized. 

Keeping these results in mind, we analyze the onset 
of PS between two non-coherent neurons of the HR type 
coupled by chemical synapses. The neurons are described 
by a 4-dimensional HR model [l3| which consists of four 
coupled differential equations: Xk — ayk + bx^ — C2:| — 

dZk +Ik+ gsyn Y^j IkjlsyniXj), ilk ^ e - yk + fxl - gWk, 

Zk = fi{-Zk + R{xk + H)),Wk = iy{-kwk + r{yk + I)), Xk 
represents the membrane potential, yk is associated with 
fast current dynamics, and {zk,Wk) are associated with 
slow currents, Igyn is the synaptic input, and jkj is the 
connectivity matrix: jkj = 1 if neuron j is connected to 
neuron fc, and "fkj = 0, otherwise, with j ^ k. 
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We set the parameters of the model in order to obtain a 
spiking/bursting behavior 13[. Then, we couple the neu- 
rons by means of chemical synapses. The current Isyn in- 
jected in the postsynaptic cell is given by Isyn{xj) — 
Sj[xj-Vrev],T^ = s^-'s^jixi) ' where Vrev is the synap- 
tic potential, and r is the timescale governing receptor 
binding. Sac is given by: SooiV) = tanh[{V-Vth)/Vsiope] 
if > Vth and otherwise 
The synapse parameters are chosen in order to have 




FIG. 3: Onset of PS in two HR neurons coupled via inhibitory 
synapses. In Figs, (a-b), we plot the attractor projection 
(xk,yk) in gray, and the set Vk in black, constructed by ob- 
serving the subsystem Efc whenever Ej crosses the threshold 
Xj = —1.3 represented by the dashed line in (c). In (a) The 
set "Dk spreads over the attractor which shows that there is 
no PS. In (b) the set T>k is localized which shows the pres- 
ence of PS. In (c) we present the time series of the membrane 
potential of (b). The threshold (dashed line) can mislead the 
burst occurrence (see the box) leading to the wrong statement 
that there is no PS. For (a) the parameters are Ik = 3.12 and 
gsyn = 0.75, for (b) and (c) Ik = 3.12 and Qsyn = 0.85. 



an inhibitory effect, so, we set: Vth = -0.80, Vsiope = 
imyrev = -1.58, and So > 1. 

Now, the time series of events rj is the time at which 
the ith crossing of membrane potential Xj reaches a 
threshold, namely x = —1.3. We fix Ik = 3.1200 and 



3.1205, then for 



0.75 the set Vk spreads over 



the attractor [Fig. [3l[a)]; there is no PS. As we increase 
gsyn, the coupled neurons undergo a transition to PS, i.e. 
the set Vk is localized. Fig. ^h). The neurons are highly 
non-coherent, due to the existence of two time-scales, and 
the inhibitory synapse which causes one neuron that is 
in a spiking behavior to inhibit the other neuron, which 
hyperpolarizes, but still tries to spike. This competition 
generates even more non-coherence in the phase space. 



As a consequence, it is rather unclear how one can cal- 
culate the phase for such dynamics. What has been cur- 
rently done is to estimate the phase of the chaotic neuron 
by assuming that in every crossing in a given direction 
of the membrane potential with a threshold, the phase 
increases 27r 0] . The main problem with this approach is 
that the phase is threshold dependent, so, it can lead to 
the false statement that PS does not exist. We illustrate 
this problem in Fig[3ljc), for the same parameter as in 
Fig[3jb); there is, indeed, PS. For a threshold x = —1.3 
(dashed line), one burst is missed, what makes the phase 
difference to be no longer bounded as the time goes to 
infinity, leading to the wrong statement that there is no 
PS. Our approach, on the other hand, is not event de- 
pendent. Indeed, as we showed, a localized set Vk exists 
for this threshold. 

Next, we analyze a network of 16 non-identical HR 
neurons, connected with excitatory chemical synapses. 
In order to simulate a mismatch in the intrinsic current, 
we set li = 3.12 -I- rji, where rji are uniformly distributed 
within the interval [—0.05,0.05]. To simulate the excita- 
tory synapses, we use the same Isyn, but changing the 
value of Vrev ■ Note that if Vrev > Xi (t) the neuron presy- 
naptic always injects a positive current in the postsy- 
naptic neuron. In the following, we set Vrev — 2. Our 
network is a homogeneous random network, i.e. all neu- 
rons receive the same number k of connections, namely 
k = 4. We constrain gsyn to be equal to all neurons. We 
identify the amount of phase synchronous neurons by an- 
alyzing whether the sets Vj are localized, occupying no 
more than 80% of the attractor of S_,- . 

The onset of PS in the whole network takes place 
gtyn ~ 0.47, so all neurons become phase synchro- 
nized. As the synapse strength crosses another thresh- 
old, gsyn ~ 0.52, the neurons undergo a transition to the 
rest state, and they no longer present oscillatory behav- 
ior. Clusters of PS appear even for gsyn far smaller than 
gsyn- In fact, right at gsyn = 0.04, some clusters of neu- 
rons exhibit PS among themselves. These clusters seem 
to be robust under small perturbations. Clusters of PS 
inside the network may offer a suitable environment for 
information exchanging. Each one can be regarded as 
a channel of communication, since they possess different 
frequencies, and therefore each channel of communication 
operates in different bandwidths. This scenario of cluster 
formation is neither restricted to this HR model nor to 
the synapse model. It can also be found in square-wave 
and parabolic bursters. 

In conclusion, we have proposed an extension of the 
stroboscopic map, as a general way to detect PS in 
coherent/non-coherent oscillators. The idea consists in 
constraining the observation of the trajectory of an oscil- 
lator at the times in which an event occurs in the other 
oscillator. We have shown that if PS is present, the maps 
of the attractor appear as a localized set in the phase- 
space, and vice-versa. The ideas herein provide a reliable 
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and easy way of detecting PS, without having to exphc- 
itly measure the phase. This method can be apphed in 
experiments in real-time and networks. 
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